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Abstract. We consider a Bose-Hubbard trimer, i.e. an ultracold Bose gas populating 
three quantum states. The latter can be either different sites of a triple-well potential 
or three internal states of the atoms. The bosons can tunnel between different 
states with variable tunnelling strength between two of them. This will allow us 
to study; i) different geometrical configurations, i.e. from a closed triangle to three 
aligned wells and ii) a triangular configuration with a 7r-phase, i.e. by setting one 
of the tunnellings negative. By solving the corresponding three-site Bose-Hubbard 
Hamiltonian we obtain the ground state of the system as a function of the trap 
topology. We characterise the different ground states by means of the coherence and 
entanglement properties. For small repulsive interactions, fragmented condensates are 
found for the 7r-phase case. These are found to be robust against small variations 
of the tunnelling in the small interaction regime. A low-energy effective many-body 
Hamiltonian restricted to the degenerate manifold provides a compelling description 
of the TT-phase degeneration and explains the low-energy spectrum as excitations of 
discrete semifluxon states. 
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1. Introduction 

It is well known that bosons at sufficiently low temperatures tend to form Bose-Einstein 
condensates (BECs), which essentially consist on the macroscopic population of a 
single-particle state [T]. In absence of interactions, the macroscopically occupied state 
is the lowest energy state of the single-particle Hamiltonian. When the interatomic 
interactions are taken into account, and for sufficiently large number of atoms, the 
main effect is a broadening of the single-particle state, which can be accounted for by 
a mean-held Gross-Pitaevskii description. In the Onsager-Penrose picture, in a BEG 
there is only one eigenvalue of the one-body density matrix which is of the order of the 
total number of particles. This is the largest eigenvalue and is termed the “condensed 
fraction”. 

In contrast, an interesting scenario appears when, in absence of interactions, the 
lowest single-particle states are degenerate. In this case, the many-body ground state 
may get fragmented [2] , as naively the atoms have no reason to condense in only one of 
the degenerate single-particle states. This implies that a hnite number of eigenvalues of 
the single-particle density matrix are of order of the total number of atoms. 

In Ref. [2] the authors describe three notable physical examples which produce 
fragmented condensates, e.g. the highly correlated regime in ultracold gases subjected 
to synthetic gauge helds, the two-site Hubbard model and the ground state of a spinor 
condensate in absence of quadractic Zeeman terms. In the hrst two cases, interactions 
need to dominate over tunnelling terms in order to get fragmentation. For instance, 
for atoms in the double-well it is in the strongly repulsive regime that the condensate 
fragments in two parts, with half the atoms populating each well. By increasing the 
potential barrier between the wells, the system enters into the Fock regime (interaction 
energy dominating over the tunnelling) and fragments in two BEGs of equal number 
of particles. This has been realized experimentally PHEIE]. Notably, the quantum 
many-body correlations present in fragmented states can End applications in the held of 
quantum metrology to improve precision measurements [S], and hold promise of being 
useful in near future technological applications [ 7 ]. 

It is desirable to pin down quantum many-body systems which feature 
fragmentation even at the single-particle level without explicit spatial separation. In 
this article we consider a minimal system which fulhlls this, and which therefore has 
fragmented ground states both in absence of atom-atom interactions or for small ones. 
We consider N identical bosons populating three single-particle states. The bosons are 
assumed to be able to tunnel between the different states. There are two options that 
are available with current techniques. The first one would be to trap the atoms in a 
triple-well potential, with fully connected sites as in ilS], or aligned |T0]. In this case, 
the three quantum states are the three eigenstates of the single-particle Hamiltonian, 
which are thus spatially localised. A second option is to consider three different internal 
states of the boson as single-particle states, with the whole cloud being trapped on the 
same harmonic potential. Transitions among the three internal states can be induced by 
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Figure 1. Schematic depiction of the system under study. We consider a gas of 
bosons which can populate three different modes, depicted as green balls. Bosons are 
allowed to tunnel between the modes. The tunnelling between modes 1 and 3 is taken 
to be tunable through the parameter 7. For 7 = 0 the three sites are aligned and the 
geometry is essentially one dimensional. For 7 > 0 we have a triangular configuration, 
which for 7 = 1 is equilateral. In the limit of 7 >> 1 the system is similar to a 
two-mode one. For 7 < 0 we have 7r-phase tunnelling between sites 1 and 3. 


means of Rabi coupling. In this case, the three quantum states are not spatially localised. 
These two cases may be referred to as external or internal, three-mode systems. The 
properties of such triple-well potentials have been studied previously exploring the effect 
of dipolar interaction in the system dn [HI El El eg, the many-body properties of the 
system [TOl ITH] or the melting of vortex states [9] . 

In the previous studies, the coupling between the different modes is provided by 
the quantum tunnelling between the spatially localised modes. This has hindered the 
exploration of the regime we discuss in the present article. Namely, we will consider 
triangular setups in which one of the tunnelling terms can be taken negative. This means 
that a particle tunnelling from one site to the next one acquires a phase of tt. Thus, we 
consider cases when a particle acquires either 0 or vr phase when tunnelling around the 
triangle in absence of interactions. The key idea is to consider systems where tunnelling 
can be detuned ini and particularly profit from the recent advances in producing phase 
dependent tunnelling terms. For external modal conhgurations, an external shaking of 
the system along one direction effectively results in a dressed tunnelling term whose sign 
can be switched from positive (standard) values to negative (7r-phase tunnelling) [TS] , 
More recently, a deep laser dip in the centre of a junction has been proposed in Ref. [19] 
to engineer vr-phase tunnelling. In the internal case a phase-dependent tunnelling can 
be obtained as in Ref. 1201 . 

This paper is organised as follows. In Sec. we present the three-mode Bose- 
Hubbard Hamiltonian, and briefly recall the many-body magnitudes that are used to 
characterise the entanglement and correlation properties in the system. In Sec. [g, we first 
discuss the non-interacting case, whose properties are crucial to understand the onset of 
fragmentation at the many-body level. Sec. is devoted to analyse the role played by 
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interactions. In Sec. |^we provide a model developed to understand the appeareance of 
fragmentation in the symmetric vr phase. Finally, in Sec.j^we summarise our conclusions 
and provide possible implications for future experiments. 

2. Three-mode Bose-Hubbard Hamiltonian 

We assume N ultracold bosons populating three quantum states. As discussed above, 
they can be different sites, e.g. an ultracold gas confined in a triple-well potential, or 
three internal states of the atoms. For the time being, we will restrict our system to the 
former case. We consider tunnelling terms between the three sites, and a tunable rate 
between two of them. This tunable rate will allow us to explore colinear configurations, 
closed ones, and also configurations with tt phase. Besides the tunnelling, the atoms 
interact via s-wave contact interactions. This interaction is assumed to be the same 
between all atoms, which is straightforward in the external case, as there is only one 
kind of atoms, but would require some tuning in the internal case depending on the 
species. The three-mode Bose-Hubbard (BH) Hamiltonian we consider is thus. 


FL 



U 


+ - 1 ), 


2=1 


( 1 ) 


where di (aj) are the bosonic annihilation (creation) operators for site i fulfilling 
canonical commutation relations, and fii = a|aj is the particle number operator on 
the Tth site. J is the tunnelling coefficient between sites 1-2, and 2-3, and U is the 
atom-atom on-site interaction that can be repulsive f/ > 0 or attractive U < 0. In 
our study, sites 1 and 3 are always equivalent with respect to site 2, and the tunnelling 
between sites 1-3 depends on the particular configuration through the parameter 7 . 

Figure shows schematically the triple-well potentials we have addressed. When 
7 = 0, no tunnelling exists between sites 1 and 3, which corresponds to an aligned triple¬ 
well conhguration. When 0 < 7 < 1, sites 1 and 3 are connected but the tunnelling 
between 1-2 (and 2-3 as well) is larger. In contrast, when 7 = 1 , the three sites are 
fully equivalent with the same tunnelling rate among them, which can be geometrically 
interpreted as arranged in an equilateral triangular potential. We will go beyond this 
symmetric conhguration, when 7 > 1 , by increasing the tunnelling rate between sites 1 
and 3 with respect to 1-2 (and 2-3), up to 7 = 2 and we will consider negative values 
of 7 besides. The last range can be engineered by lattice shaking along the direction of 
sites 1 and 3 [IB]. 

The Hamiltonian of Eq. ([^ can also be reproduced with a three-component 
(spinorial, isotopic or atomic) BEC mixture trapped on a single harmonic oscillator 
with suitable Rabi coupling between the levels (which corresponds to the tunnelling 
terms). 

Our study will concentrate mostly on the repulsive interaction case, which is more 
prone to be experimentally explored with current setups. We will also discuss briehy 
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the attractive interaction case. In the latter even small asymmetries in the external 
trapping potentials will eventually have a large impact on the properties of the ground 
state izu. Thus, we will introduce small symmetry breaking terms, compared to both 
the tunnelling and interaction, to consider situations closer to experimental ones and 
also to control numerically the degeneracies in the problem. We have added the biases 
between sites 1-3 and 2-3, eis and 623 , respectively. The considered Hamiltonian now 
reads: 

ii ' = ii + ei3(hi - hz ) 633(^2 - hg). (2) 


2.1. Many-body basis 


The system is studied through direct diagonalisation of the many-body Hamiltonian for 
a hxed number of atoms, N. A suitable many-body basis is the Fock one, which labels 
the number of atoms in each mode. 


ni,n2,n3) 


1 



vac) , 


(3) 


where |vac) stands for the vacuum, and N = ni + n 2 + n^. The elements of the Fock 
basis can be expressed as a product state |ni,n 2 ,n 3 ) = |ni) ® \n 2 ) ® 1 ^ 3 ). A general 
many-body wavefunction is thus written as 


N 

Cni,n2 1^15^25^3) 1 ( 4 ) 

ni,n2 


where Cni,n 2 is fhe corresponding amplitude of the Fock state |ni,n 2 ,n 3 ), and = 
N - {ni+ 712 ). 


2.2. Coherent states 

Coherent states are the closest analogs to classical solutions, in the same way 
wavepackets are the closest quantum analog to classical trajectories. A general coherent 
state can be constructed by assuming that all N atoms populate the same single-particle 
state i, 

h\ = cid\ + C 2 cil + 0303 . (5) 

The coherent state reads, 

I'J'coh)=^^ I vac). (6) 

Since q G C, this many-body state has 6 parameters to be determined. Properly 
normalising the single-particle wavefunction and also realising that there is always an 
arbitrary global phase, the number of free parameters can be reduced to 4. It is worth 
stressing that a coherent state as the one dehned above corresponds to a fully condensed 
atomic cloud. 
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The fragmentation properties of the ultracold atomic gas |2] can be investigated by 
means of the eigenvalues of the one-body density matrix. In our system with N bosons 
and three different modes, the one-body density matrix of a many-body state |\h) is a 
3x3 matrix whose elements are, 


( 1 ) 

u 


1 

N 


(T|ala,|T), 


(7) 


with i,j = 1,2,3. Since |T) is normalised, it follows that Tr= 1. 

It is interesting to calculate its eigenvectors (natural orbitals), |'0i), and eigenvalues. 
Pi, with Pi > P 2 > Ps > 0, that satisfy Pi -l- P 2 + Ps = 1- Each eigenvalue of the one- 
body density matrix gives the relative occupation number of the corresponding natural 
orbital: pi = Ni/N. In a singly condensed system, there is only one large eigenvalue 
that corresponds to the condensed fraction of the single-particle state, \'ipi). Pi ~ 1 and 
A^i ~ 0{N), with all the other eigenvalues pj {j ^ 1) being small ~ 0{1/N). Instead a 
fragmented system has more than one large eigenvalue, Ni ~ 0{N), with z = 1,..., s, 
and the rest of eigenvalues pj {j > s) are small ~ 0{1/N). In this situation the system 
is not fully condensed, but fragmented, quantum correlations become important and 
the mean-field approximation fails to describe the system. 


2 . 4 . Entanglement properties: entanglement entropy and Schmidt gap 


Correlations between different subsystems of a many-body quantum system can be 
quantified performing different bipartite splittings. That is, considering the system as 
made of two subsystems, tracing out one of the parts, and studying the von Neumann 
entropy and the entanglement spectrum [22] of the resulting subsystem. 

In our case, we consider different spatial partitions of the three-well conhguration, 
e.g. (1,23), (2,13), as in Ref. (H]. From the density matrix of the full system, p, 
correlations between mode i and the rest can be determined by first taking the partial 
trace of p over the Fock-state basis of the other modes. This yields the reduced density 
matrix on subsystem i, Pi, that describes the state of this subsystem. For instance, in 
our system, by tracing out sites 2 and 3, a bipartite splitting of the three-mode system 
is obtained (1, 23), and the reduced density matrix on site 1, pi = Tr 23 p, is found to be 


diagonal in the single mode space of N particles (see Eq. (A.2) in Appendix Appendix 

0 ’ 


N 


Pi = z2^k\k){k\, A^>0 

k=0 


( 8 ) 


where \k) are states of k particles in mode 1. Note, that the reduced density matrix 
for state 1 is in general a mixture without a well-defined number of particles. The 
set of eigenvalues {A^} is called entanglement, or Schmidt spectrum [2211231 [21|, and 
the eigenvalues are the Schmidt coefficients. The Schmidt coefficient is in this case 
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directly the probability of finding k particles in site 1 withont measnring the nnmber 
of atoms in sites 2 and 3. The Schmidt spectrnm fnlfills = 1, and contains 

information abont the correlations and the entanglement properties of the state in 
snbsystem 1 with respect to the rest of the system. It is worth recalling that a many- 


body state is entangled when it cannot be written as a prodnct state. In the case of 


spatially separated modes, the entanglement we are discnssing is spatial. 

A measnre of the entanglement between the two snbsystems is already provided by 
the single-site von Nenmann entropy, which can be compnted as Si = —Tr(pi logpi). 
Noting that in our case the density matrix, pi, is already diagonal, see Eq. (|^, the 
entropy can be evaluated from the Schmidt coefficients as S'! = — Aj log A^l^. 

In the three-mode system, a signature of the entanglement on the Schmidt spectrum 
is the following one: if site 1 is not entangled with sites 2 and 3, it should be pure after 
tracing those sites out, and then the entanglement spectrum would be A} = 1 and 
Aj = 0 for i = 2,3,.... This actually implies a zero of the corresponding von Neumann 
entropy. A remarkable magnitude defined from the set of A’s is the so-called Schmidt gap, 
defined as the difference between the two largest and more relevant Schmidt coefficients 
in the entanglement spectrum of subsystem i, AA* = A| — A 2 |23]- In the case of no 
entanglement between the subsystems, the Schmidt gap takes its maximum value of 1. 
A vanishing of the Schmidt gap marks large entanglement between the subsystems. As 
has been recently pointed out in Ref. na for dipolar bosons in triple-well potentials, 
the Schmidt gap is a good tool to distinguish between phase transitions and crossovers. 

Note that no relation exists between fragmentation and entanglement, in the sense 
that a system can be entangled and fragmented independently. For instance, a system 
where all the bosons occupy the same spatial mode is nor fragmented neither spatially 
entangled. A bosonic Josephson junction in the strong repulsive interaction (Fock) 
regime is a clear example of a non entangled but fragmented state: |iV/2,A/2). In 
contrast, in the non-interacting regime the bosonic Josephson junction is in a fully 
condensed state. 



(9) 


which has large entanglement between the two sites as seen by the Schmidt coefficients 
which are Xk = 2“^('^). Finally, cat states, |T) = (|iV, 0) -|- |0 ,A))/a/ 2, are a well 
known example of fragmented and entangled many-body systems. 

3. Non-interacting case 

In this section we solve the non-interacting case for any value of the tunable tunnelling 
link 7 J. In our calculation we will keep the tunnelling parameter fixed J/h = 1 Hz, 
since it essentially sets the overall energy scale. 

The triple-well configuration we have chosen is intended to remark first, the role 
played by the topology of the configuration, going from a disconnected chain (7 = 0 ) to 
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Figure 2. Single-particle spectrum of the triple-well system as a function of 7 . In 
black, red and green solid lines we depict the eigenvalues Ei , E 2 and E ^, from Eqs. (10 1 . 
The corresponding eigenstates, Eq. (11), are written explicitly for several relevant 
values of 7 . The dashed lines correspond to the spectrum computed with a bias term 
£13 = £23 = <^/4- The short notation for states (a, b, c) in the figure corresponds to the 
states a|l, 0,0) -I- 6|0,1, 0) -I- c|0,0,1) used in the main text. 


an essentially donble-well system at 7 —)■ cx), throngh a connected eqnilateral triangle at 
7 = 1. It is worth emphasising that we do so by varying the tnnnelling strength between 
sites 1 and 3 and not by, for instance, adding sizeable bias terms to the Hamiltonian, 
which wonld indeed also break the symmetry between the three sites. The second 
important point, is that we consider vr-phase tnnnelling between the two wing sites 
(7 < 0), which will indeed have dramatic conseqnences on the many-body properties of 
the system. 

3.1. Unbiased non-interacting case 

The non-interacting case can be solved readily. The eigenvalnes of the Hamiltonian (|^ 
with U = 0 are. 


El = 7J 

E 2 = ^ (-7 - 7^8 + 72^ 

Es = ^ (-7 + ^8 + 7^) , 

and the corresponding nnnormalized eigenvectors are: 


( 10 ) 
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|V> 1 > = ^(| 1 . 0 , 0 >-| 0 , 0 , 1 )) 

Ife) - 11.0,0) + ^ ~ 10.1,0) + |0,0,1) 

v 8 + 7 ^ + 37 

I Va) ^ | 1 . 0 , 0 ) + — |0,1,0) + |0,0.1). (11) 

V 8 + 7 ^ - 37 

The single-particle spectrum, see Fig. has some interesting features. First, 
there is one eigenvector that is independent of 7 , \4>i)- This state does not involve 
site 2, and its energy is proportional to the coupling between sites 1 and 3. For 
7 = 0 , which corresponds to an aligned conhguration, the ground state is 1 - 02 ) 7=0 = 
1/2 (| 1 , 0 , 0 ) -|- a/2 |0, 1 , 0 ) -|- | 0 , 0 , 1 )), which has an excess of particles in site 2 . 

As 7 is increased, the equilateral triangular conhguration is reached for 7 = 1 . 
This case has been studied thoroughly in Refs. PEg. The ground state is |' 02 ) 7 =i = 
(l/-\/3)(|l, 0 , 0 )-|-| 0 ,l, 0 )-|-| 0 , 0 ,l)) and there are two excited degenerate single-particle 
states that correspond to vortex states, with clockwise and counter-clockwise rotation: 

10.1) = ^ (|1,0,0) -F |0,1, 0) + |0, 0,1)) 

10 . 2 ) = ^(| 1 , 0 , 0 )+ e-*^| 0 ,l, 0 )+e-*^| 0 , 0 ,l)) . ( 12 ) 

Notice that these two states are linear combination of 

iV’i) = ^ (e"*^ 10 ^ 1 )+ 10 ^, 2 )) 

103)7=1 = 0’ 0) - 2 | 0 , 1 , 0 ) + 10 , 0 , 1 )) 

= ^ (e*^ 10^1) + |0„2)) . (13) 

Further increasing 7 , sites 1 and 3 get further connected and the physics decouples 
them from site 2. The ground state for 7 —)■ cx) is 102 ) 7^.00 = l/v^ (|1, 0, 0) -|- |0, 0,1)), 
the hrst excited state is | 03 ).^^oo = |0,1,0) and the second excited state is | 0 i) 7 ^oo = 
1 / 72 ( 11 , 0 , 0 )-| 0 , 0 , 1 )). 

The situation for 7 < 0 is very different and actually richer at the ground- 
state level. For 7 = —1 we have a crossing in the single-particle spectrum, which 
therefore should have important consequences at the many-body level. At 7 = 
—1 the ground state is two-fold degenerate between states | 0 i) and | 02 ) 7 =_i = 
(|1, 0, 0)-|-2|0,1, 0)-|-|0, 0,1)) /a/6. The hrst excited state is = (|1,0,0) — 

|0,1, 0) -|- |0, 0, 1))/a/3. When varying 7 from 7 < —1 to 7 > —1, the system undergoes 
a quantum phase transition in this non-interacting case. A possible order parameter 
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would be the population of the second site, which goes from (72,2) = 0 to (722) > 0 as we 
cross the 7 = — 1 point. 

From Fig. one can see that the main effect of varying 7 is a large avoided crossing 
of the asymptotic states (|1,0,0) + |0,0, 1))/a/2 and |0,1,0). When 7 is positive the 
ground state is (|1,0,0) + |0,0, l))/\/2, whereas when 7 is negative the lowest state is 
the antisymmetric state, l-^i) = (|1,0,0) — |0,0, 1))/a/2. The latter remains uncoupled 
for all values of 7, with eigenenergy E = J7. This value is easy to obtain, since in the 
Fock basis: [a\a 3 + a\a\\\'tp\) = Therefore, {'ipi \ — + a3ai)|V’i) = J'l- 

3.2. Explicit symmetry breaking, effect of the bias 

The explicit symmetry breaking induced by bias terms in the Hamiltonian (|^ breaks 
the degeneracies present in the single-particle spectrum. As seen in Fig. the crossings 
at 7 = — 1 and 7 = 1, which occur in the non-interacting system for the ground state 
and excited states, respectively, become now avoided crossings. For the case we are 
interested in, when the bias is much smaller than the tunnelling, we can obtain the 
states dressed by the bias at the degeneracy points by means of hrst order perturbation 
theory. In the case of 613 > 0 and 623 = 0, the ground-state manifold at 7 = — 1 splits, 
AE = ei3, and the corresponding dressed states are. 


1^1) = ^ (1^2) - iV’i)) = 


V12 


(1 - \/3)|l, 0,0) + 2|0, 1, 0) + (1 + v^)|0, 0, 1) 


1^2) = ^ (IV'i) + IV'2)) 


^ L 


(1 + 73 ) 11 , 0 , 0 ) + 2 | 0 , 1 , 0 ) + (1 - 73 ) 10 , 0 , 1 ) 


(14) 


The hrst one is the new ground state of the system, which has a slightly larger population 
of particles in site 3. This result is reasonable, since the bias we have considered, 613 > 0, 
promotes site 3. 


4. Quantum many-body properties of the system 

We consider now the effect of repulsive interactions between atoms. We calculate the 
ground state, by exact diagonalisation of the Hamiltonian ([^, for different values of 
the tunnelling rate 7. In our numerics we will consider up to = 48 particles. The 
ground state of the system is characterised by means of, a) coherence properties and 
fragmentation, b) analysis of the energy spectrum, and c) entanglement spectrum and 
entanglement entropy. 
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Figure 3. Phase diagram in the ( 7 , NU/ J) plane obtained by exact diagonalisation for 
N = 48. The colour gives the average occupation in the ground state, (4'Gsl^il^Gs)/-^) 
of site 1 (a) and 2 (b). 


Analytic results in the I7I S> 1 limit 

The structure of the single-particle spectrum, see Fig. allows one to build a simple 
model for I7I 1. This simple model will serve as guidance to understand many of the 
properties of the many-body ground state which will be discussed in the forthcoming 
sections. 

We can distinguish two distinct regimes, a) 7 S> 1 and b) —7 S> 1. In both cases, 
for sufficiently large values of I7I the single-particle spectrum approaches Eq = — |7|J, 
El = 0, and E 2 = \'y\J- In the (a) case, the non-interacting single-particle states are 
iV^o) = 1/72(11,0,0) + |0,0,1)), iv^i) = |0,1,0) and iV-a) = 1/72(|1,0,0) - |0,0,1)), 
respectively. In (b) the single-particle states exchange their role, with 1-02) being 
the ground state and |'0o) the more energetic. For interaction strengths such that 
NU / J -C I7I, the effect of the interaction can be considered perturbatively. In this limit 
we can ignore particle-hole excitations to the highest single-particle state, and write 
configurations with k excited atoms as, 

l'r)=A/'(4.)''^'‘(7.)‘l''“=) (15) 

where A/" is a normalisation constant and creates a particle in the single-particle 
state |'0j). The energy for this state, up to constant terms, reads, 

AT — jb/Af — h-l-1 \ TJ 

E{N, k) = U^- ^ - 1 j + -k{k - 1) + \j\Jk . (16) 

In the non-interacting case, the energy is minimal for A) = 0, as expected. For non-zero 
interactions, this expression has a minimum for 

'N f 

3 V NU )\’ 


k = Int 


(17) 
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Figure 4. Condensed fraction pi (largest eigenvalue of the one-body density matrix) 
as a function of 7 and NU/J. Notice that in this parameter regime pi + p 2 — 1 and 
P 3 ~ 0. N = 48. 


where the Int[...] stands for the integer part of the expression. Incidentally, for large 
interactions, NU/J ':J> 1 (with \'y\J 3> U), the minimal energy is obtained for k = N/3, 
as expected for a Mott insnlator of hlling N/3. As we decrease interactions, the system 
goes step by step to A; = 0. In particnlar, many-body states with k and A; -|- 1 excited 
atoms degenerate if. 


NU _ 2N\'y\ 

~T ~ N - 3k- 1 ' 


( 18 ) 


In this limit, it is qnite reasonable to expect that as we increase the atom-atom 
interactions keeping 7 fixed, the system tends to minimise the nnmber of pairs per site, 
which is achieved by eqnipopnlating the three sites. Eq. ( |I^ predicts gronnd state 
degeneracies, i.e. energy crossings, for certain valnes of I7I and NU/J. 

This model works reasonably both for 7 < — 1 and 7^1, where already the 
strnctnre of the single-particle spectrnm starts to resemble the asymptotic one. There 
is one important difference between 7^—1 and 7^1. In the former case, the lowest 
energy single-particle state is independent of 7 (see Fig. and has no popnlation of 
site 2. This makes the model ontlined above fairly accnrate to describe the different 
transitions in the many-body gronnd state. In the latter however, the lowest single¬ 
particle state only asymptotically approaches the state (l/\/2)(|l, 0, 0) -|- |0,0,1)). In 
this case, the model only provides a qnalitative pictnre at smaller valnes of |7|. Note 
also, that for I7I 3> 1, the valne of k which minimises the energy, corresponds in this 
limit to the average popnlation of site 2 in the gronnd state, while the average popnlation 
of sites 1 and 3 is {N — k)/2. 
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Figure 5. The two largest eigenvalues of the one-body density matrix, pi (solid) > 
P 2 (dashed), are depicted as a function of 7 (left) and NU/J (right). 

4-2. Coherence and one-body density matrix 

To have an overall picture of the different regimes that we will encounter for different 
values of NU / J and 7 we hrst study the average populations of the three sites, 
(TGg|fij|TGs)/Af. Notice that in our system, without bias terms, sites 1 and 3 are 
equivalent. In Fig. |^we present results for —2 < 7 < 2 and relatively small values of 
the interaction NU/J < 10, which do not reach the Fock regime. For very small values 
of the dimensionless parameter NU/J, one does not expect substantial changes from 
the single-particle case. Indeed in the ground states for 7 > — 1 the three modes are 
substantially populated, whereas for 7 < —1 the second mode is clearly less occupied 
than the other two as reflected in the dark region on the left corner of Fig. (b). 

In the non-interacting limit, U = 0, the problem becomes a single-particle one and 
in the case of a symmetric conhguration, 7 = 1 , the many-body ground state can be 
written as [a [25], 

= —f- 

in which the average population on each site is, for symmetry reasons, N/3. 

For large repulsive atom-atom interactions, regardless of the value of 7 , the system 
will fragment in an effort to diminish the number of pairs inside each site, in analogy to, 
for instance, the double-well |2|- In the large interaction limit, that is in the Fock regime 
t/ 3> J, the ground state can be well approximated by IH/gs"'^) — 1-^/3,-^/3, iV/3), 
which is the equivalent to a Mott insulator of hlling N/3. This ground state is non¬ 
degenerate if the number of bosons is commensurate with 3. If the number of bosons is 
not proportional to 3 the ground state becomes three-fold degenerate in the strong Fock 



+ ^2 + ^3] 


Ivac), 


(19) 
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Figure 6. Depiction of the three ordered eigenvalues of the one-body density matrix as 
a function of NU/ J for 7 = —1. The figure shows the transition from bifragmentation 
in the NU/J ~ 1 regime, to trifragmentation in the NU/J —)■ 00 (Fock) regime. 

regime: |A^/3 ± 1, N/3, N/3), \N/3, N/3 ± 1, N/3), and \N/3, N/3, N/3 ± 1), where the 
plus (minus) sign refers to a single particle (hole) delocalization. 

In the Josephson regime, defined as NU/J ~ 1, the ground state of the system is 
mostly condensed for 7 < —0.5. This is reflected in the eigenvalues of the one-body 
density matrix, see Fig. |^and left panel in Fig. As already pointed out in the single¬ 
particle spectrum, a very interesting feature is readily found in the vicinity of 7 = — 1 . 
In this case, the ground state of the system is fragmented in two pieces even in the 
non-interacting case. As the interaction increases, but still in the Josephson regime, 
the system is seen to remain bifragmented. From Fig. one can see that the region 
in the 7 -space, around 7 = — 1 , which corresponds to a fragmented condensate, slightly 
increases with the interaction. 

As discussed above, in the limit of the strong Fock regime, the ground state should 
be essentially fragmented in three pieces (in our case N = 48, which is multiple of 3), 
corresponding to the ground state |A^/3, N/3, N/3). This is clearly seen in Fig. which 
depicts the behavior of the three eigenvalues of the one-body density matrix in the 
7 = — 1 conhguration for different values of the interaction. In the non-interacting case, 
the system is fragmented in two condensates (pi = P 2 = 0.5 and ps = 0), whereas as 
NU / J increases, pi = P 2 < 0.5 decrease and ps > 0 increases, fulfilling pi + P 2 + Ps = 1- 
Moreover, one can see that Pi —?• 1/3 asymptotically. Note however that the origin 
of bifragmentation is directly related to the degeneracy at the single-particle level and 
remains in the presence of tunnelling. Trifragmentation requires strong interactions 
such that essentially tunnelling plays no role and the system can be regarded as three 
independent condensates. 
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Figure 7. (a) Energy gap between the ground state and the first excited state, (b) 
Energy difference between the second and hrst excited states. N = 48. 


In the right panel of Fig. the behavior of the two largest eigenvalues of the 
one-body density matrix is shown, as a function of NU/J, for a hxed 7 conhguration. 
At 7 = —0.99 there is a sharp transition from a singly condensed (at = 0 ) to a 
bifragmented system already for very small values NU/J > 0. However, as 7 departs 
further from —1 (7 = —0.7 and —1.3) the transition between both regimes becomes 
smoother, and the system remains fully condensed for a larger range of interactions. 


4-3. Energy spectrum 


The model presented in Sec. 4.1 predicts a number of degeneracies in the many-body 


energy spectrum of the system. Indeed, Eq. (18) predicts iV/3 energy crossings for 
both 7 < —1 and 7 > 1 . As explained above, these predictions are expected to hold 
for 7 < —1 and are indeed found in the exact many-body spectrum as seen in Fig. 


Varying the value of k in Eq. (18) from 0 to iV/3 one obtains the large I 7 I behavior of 
the different lines of zero gap in Fig. (a). Similarly, albeit not shown in the hgure, for 
7 1 the corresponding gapless lines are also found in the system. 

For 7 > 0, also interesting is the closing of the energy gap between the hrst and 
second excited states for 7 = 1 (see the vertical line in the right panel of Figs. [^. It 
corresponds to the degeneracy between vortex and antivortex states studied in Ref. [9] , 
whose wavefunctions have been previously obtained in the non-interacting case, see 


Eq. (12). As the interaction is increased the degeneracy between the corresponding two 


states remains. 


4-4■ Entanglement spectrum and entanglement entropy 

The many-body ground state has been found to be fragmented in the vicinity of 7 = — 1 
and mostly condensed otherwise, for NU/J ~ 1. Besides fragmentation, the three-site 
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Figure 8 . Schmidt gap as a function of ( 7 , NU/J) for N = 48. (a) Between site 1 
and the subsystem formed by sites 2 and 3. (b) Between site 2 and the subsystem 
formed by sites 1 and 3. 


configuration considered allows one to study the onset of entanglement and correlations 
among the three different sites depending on the specihc values of the parameters. To 
characterise the spatial entanglement properties we will use the Schmidt gap and the 


entanglement von Neumann entropy defined in Sect. 2.4 


Due to the structure of our system, with three sites where two of them, 1 and 
3, are essentially equivalent, one can consider two bipartite splittings. The hrst one 
corresponds to subsystem 1 where sites 2 and 3 have been traced off (1,23). And 
second, subsystem 2 where subsystems 1 and 3 have been traced off (2,13). 

The Schmidt gap corresponding to both bipartite splittings in the {'y,NU/J) 
diagram is presented in Fig. The hrst notable feature is that in both bipartite 
splittings we observe two fan-like radial structures, with straight lines converging to 
(—1,0) and (l,c) (with c a constant to be determined later). This structure, similar 
to the one discussed in Ref. mi, represents the crossings of the hrst two values of 
the entanglement spectrum. To better understand the structure, in Fig. [^we plot the 
entanglement spectrum [N +1 coefficients) for both bipartite splittings for a hxed value 
of NU/J = 2. For 7 < —1, each zero of the Schmidt gap shown in Fig. implies a 
variation of one unit of the most probable population of the untraced mode. This can 


be understood from the expression in Eq. (A.2), where the reduced density matrix is 


shown to be diagonal in the Fock basis of the untraced mode. 


The simple large I 7 I model explained in section 4.1 also explains qualitatively the 


observed behaviour. The degeneracies in the energy spectrum obtained there correspond 
to many-body ground states in which the most likely value of the population in site 2 
goes from 0 to N/?> and the corresponding one of sites 1 and 3 from N/2 to A^/3. We 
have considered N = 48, see Fig. and thus we find 16 crossings in the left part of the 
hgure, with the most likely population of site 1 going from N/2 = 24 in the leftmost 

























Fragmented condensation in Bose-Hubbard trimers with tunable tunnelling 


17 




y 

Figure 9. Entanglement spectrum of the reduced density matrix of mode 1 and mode 
2, upper and lower panel, respectively. The spectrum is coloured according to the 
corresponding most probable occupation of the mode singlelling out some relevant 
cases, e.g. empty mode fc = 0, half occupation k = N/2, etc. Both are computed for 
NU/J = 2 and TV = 48. 


case to N/3 = 16 at 7 = —1. As explained in the previous section, the perturbative 
model captures the physics also for 7 > 1 . As we can see in the hgure, at 7 = 1 the most 
likely population of site 1 is again N/3 = 16 and it keeps increasing as 7 is increased. 

These crossings described in the entanglement spectrum, clearly visible in Fig. 
are one of the main signatures of a crossover between two phases, whereas in a quantum 
phase transition all the Schmidt coefficients degenerate to the same value in the 
thermodynamic limit {N —)■ cx)) [23]. In our case, it is clear that at (7 = —1, NU/J = 0) 
all the Schmidt coefficients degenerate even for finite N. This is due to the fact that we 
already have a single-particle degeneracy. 

Certain limiting cases are easily interpreted. For instance, for 7 <C —1, mode 
2 is essentially unpopulated and decoupled from modes 1 and 3. This reflects in a 
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Figure 10. Entropy as a function of ( 7 , NU/J) for N = 48. (a) Bipartite von 
Neumann entropy between site 1 and the subsystem formed by sites 2 and 3. (b) 
Bipartite von Neumann entropy between site 2 and the subsystem formed by sites 1 
and 3. 


singly populated entanglement spectrum for NU/ J = 2 in Fig. (lower panel). In this 
regime, the system was shown to be condensed in Fig. on the single-particle state 
(11, 0, 0) -b |0, 0,1))/-\/2, which spatially entangles modes 1 and 3, as seen in Fig. (upper 
panel). In the 7 3> 1 limit the situation is exactly the same, but as explained above, this 
limit is achieved in practice for much larger values of I 7 I than in the 7 < 0 case. In the 
vicinity of 7 = 0 the three single-particle states quasidegenerate. This makes that for 
relatively low interactions, as in Fig. the many-body ground state starts populating 
the Fock states around N/?> approaching the Mott insulator phase. 

These features, described on the full entanglement spectrum, reflect directly on 
the corresponding von Neumann entropies, see Fig. For instance, the fact that the 
system approaches the Mott regime for relatively low values of the interaction in the 
vicinity of 7 = 0 reflects almost zero value of the von Neumann entropy for mode 1, 
which starts to decouple from the other modes. For 7 < — 1 and small values of the 
interaction, the system empties mode 2 and decouples it from the other two modes, see 
Fig. 10 (right panel). 


4-5. Quantum phase transition for attractive interactions 

Similarly to the two-well system, a quantum phase transition can be expected for 7 = 1 
for attractive interactions [291130]. In this case, the three sites are completely equivalent. 
Increasing the attractive atom-atom interactions, the system will minimise energy by 
clustering the atoms in a single site. In absence of any spatial bias, the ground state of 
the system will approach the Schrodinger cat-like state, 

|>P„) = 4(N, 0 ,0) + |0, N, 0) + |0,0, ]V». 


( 20 ) 
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Figure 11. Schmidt gap as a function of ( 7 , NU/J) for N = 48. (a) Between site 
1 and the subsystem formed by sites 2 and 3. (b) Between site 2 and the subsystem 
formed by sites 1 and 3. This figure is similar to Fig. but we extend to attractive 
interactions to show the phase transition described in the text. 


Analogously to the two-site case, this many-body state is not gapped and is 
quasidegenerate with its hrst two excitations. In the thermodynamic limit, the transition 
between the non-interacting state, Eq. (19), and the cat-like state Eq. (20) goes through 
a transition point at a hnite value of NU / J. This transition reflects in the behaviour 
of the Schmidt gap of the system. As seen in Fig. several zero-Schmidt-gap straight 
lines tend to converge on a point in the attractive interaction regime on the 7 = 1 line. 
In Fig. 11 we extend the range of parameters to the attractive region and certainly the 
lines seem to converge at 7 = 1 and —9/2 < NU/J < —4, which is where the authors 
of Ref. [13] predicted the existence of a quantum phase transition. Notably, in contrast 


with the former phase transition described at 7 = — 1 , this phase transition is only 
present in the thermodynamic limit. 


5. Detailed analysis of the 7 = — 1 case 

As already shown in Figures and the 7 = — 1 ground state is found to be 
bifragmented, with pi = P 2 = 1/2, in the non-interacting limit {NU/J = 0). To a good 
approximation it remains bifragmented for hnite but small {NU/J ~ 1) interactions. 
As NU/ J increases further the system approaches the Fock regime and the ground state 
tends to the well known trifragmented conhguration with pi = P 2 = Ps = 1/3. We will 
here instead focus on the description of the special bifragmented states obtained when 
NU / J is non-vanishing but small. 

At 7 = — 1 the single particle spectrum, see Fig. has two degenerate eigenvalues. 
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Es = Ea = —J and E^ = 2 J. A convenient choice for their eigenvectors is, 

1 ^) =^(| 1 , 0 , 0 )+ 2 | 0 , 1 , 0 )+ 10 , 0 , 1 )) 

|A) = ^(|1,0,0)-|0,0,1)) 

|e) =^(|1,0,0)-|0,1,0)+ 10,0,1)), (21) 

so that, under the exchange of sites 1 O 3, \S) is symmetric and |A) is antisymmetric. 
|e) corresponds to the “excited” single-particle conhguration. Eg — Es = 3J. To these 
states we associate creation operators 

= —1= + 2(32^ + 


ye 
V2 




e* = ^ (aj - aj + a: 


( 22 ) 


and their hermitian conjugates, S, A and e, act as annihilation operators. It is immediate 
that 

+ + = 1 , [y+ = l , [e,et] = l. (23) 

Formally, it is a matter of convenience to use modes 1, 2, and 3 or S', A and e to 
construct the single-particle basis. 

To understand the onset of bifragmentation we will now develop an approximate model 
in the following way: we will only use the subspace generated by the and modes 
to describe the lowest lying many body states, in particular the ground states. To do 
so, inverting Eq. (22) one hnds aj in terms of A^, S'^ and e^ but since the latter do not 
operate in the chosen subspace, in our approximation, 

at = ^At + ^,Ft 
y/2 ye 

dl = \[^S^ 


i\ = - ^At + ^At . 


V2 ye 


(24) 


Within that subspace the interaction Hamiltonian reads, 

3 

^ t \ ^ 

a, 


flv = E+)+ = t(i*)"i" + 

i=l 

2 1 


+ ^ (^1)2^2j ^ ^25) 

where for simplicity we have omitted the U/2 factor and the terms linear in N. 
Introducing two new modes. 


A _ 


^hvl 
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Figure 12. Overlap between the exact two-fold quasidegenerate groundstate manifold, 
!$_,) and the two k = 0 states in Eq. (30). The overlap is computed as |detO|, with 
Oij = i = +,— and j = 1, 2. Note this overlap is 1 if the two sets of vectors 

span the same manifold. 


a 


L = - iA<) 


the Hu Hamiltonian becomes, 
1 
3 

- at 


Hu — - ( -|- Nh,o. + 4:NhviNhv2 — Nhvi — Ni 


hv2 


(26) 


(27) 


hv2 

where = d\^-^^dhvi and Nhy 2 = (ihy 2 ^hv 2 - From Eq. (27) it is clear that the Fock 
states built from the new modes 

1 


\Nhvl, Nhv2) = 


VNhvl^-Nhv2^- 


(“Li)"*”' (“’m)"*'"!™':) 


(28) 


are the eigenvectors of Hu, with an obvious expression for the eigenvalues. Since 
Nhvi + Nhv 2 = N, it will be more clarifying to change notations to Nhvi = k and 
Nhv2 = N — k. Then the spectrum becomes 

Eu{k) = \{k^ + {N-kf + Ak{N-k)-N) , (29) 

o 

which is degenerate in pairs, {k,N — k), except the topmost energy when N is odd. In 
particular, the ground state is two-fold degenerate. 

This approximation turns out to be very accurate for small interactions. The effect 
of including the |e) manifold, breaks the degeneracy but does not promote any of the 
states. The lowest many body states are then well approximated by. 
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Figure 13. Quantum phase in each site, color value, for the fluxon Eq. (121 and 
the semifluxon Eq. (311 states. The fluxon is the usual vortex state, with a full 27r 
variation of the phase around the triangle. In the semifluxon, the phase varies by tt as 
we go along the triangle. 


with k = Q,...,N 


(30) 


where the ± sign labels the two states which would be degenerate in energy in absence of 
coupling with the |e)’s. The are obviously bifragmented, and have pi = P 2 = 1/2. 
Furthermore, the ground state is a cat-state with all atoms in one or the other of the 
modes in Eq. (26). These can be considered creation operators of discrete semifluxon 
states [32] because, combining Eqs. (22) and Eq. (26), they can be written as 

1 


at - 

^hvl ~ 


‘'hv2 


Vs 

1 


+ dl 


+ e 


-W3^t 


= —jzz ( e - 1 - 02 - 1 - e' 


t 


T/S-t 


(31) 


A semifluxon is a quantum state which carries half the quantum flux of the vortex states 
of Eq. (12). In our case, we have a discrete version, as going around the triangle the 
quantum phase grows from 0 to tt, with the phase jump of vr imposed by 7 = —1. In 
Fig. 13 we depict the phase structure of the discrete semifluxon state, compared to the 
usual vortex, Eq. (12) (referred to as “fluxon”). 

In the full three-mode space, as said above, the degeneracy of each doublet 
splits and therefore the ground state should be a NOON state of both semifluxon states. 
In Fig. I^we depict the overlap between the exact ground state and hrst excited state 
and the corresponding /c = 0, ± states in Eq. (30). The model is seen to be fairly 
accurate up to NU/J ~ 10. In an experimental measurement, one should thus hnd a 
fractional flux in one or the other direction. 
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6 . Summary and conclusions 

A simple configuration consisting of three sites with a single tunable link has been shown 
to exhibit a large variety of quantum many-body properties. An important novelty of 
the proposed scheme is that we have considered cases in which the tunnelling rate in one 
of the links is either zero-phase or vr-phase tunnelling. In both cases, we have performed 
diagonalisations for hnite number of particles, up to = 48, and have examined the 
many-body properties of the systems as a function of the atom-atom interaction and 
the tunable tunnelling rates. 

By varying the tunnelling in one link, by means of the parameter 7 , which is 
within reach experimentally as explained in the introduction, one explores conhgurations 
ranging from the colinear one 7 = 0 , the fully symmetric one 7 = 1 , to the symmetric 
TT-phase one, 7 = —1. In the hrst two cases, the many-body ground state for small 
interactions NU / J ~ 1 is highly condensed, with one of the three eigenvalues of the one- 
body density matrix clearly scaling with the total number of particles. As we approach 
the 7 = — 1 point the system departs from condensation and becomes bifragmented for 
small interactions. This is partly a consequence of the degeneracy in the single-particle 
spectrum for 7 = — 1 . 

Two quantum phase transitions are present in this system. The hrst one takes 
place in the symmetric conhguration 7 = 1 , but only for attractive interactions, 
which therefore makes it difficult to explore experimentally. This phase transition 
is reminiscent of the one present in other few-site models, e.g. bosonic Josephson 
junctions, and marks the transition to Schrodinger cat-like ground states in the 
spectrum [221 1321 El]- Interestingly, for 7 = —1 we hnd a second quantum phase 
transition, which takes place in absence of interactions and which has clear consequences 
for small repulsive atom-atom interactions. This phase transition has been characterised 
by the behaviour of the entanglement spectrum for the two possible independent 
bipartite splittings as we approach the transition point [23]. At the transition the 
entanglement spectrum degenerates and the corresponding von Neumann entropies 
exhibit a maximum. Interactions do not wash out the main features of this transition, 
which becomes essentially a crossover, which has clear consequences at hnite interaction, 
such as the bifragmentation discussed above. 

Finally, for the 7 = — 1 case and small interactions, we have been able to derive 
an approximate many-body Hamiltonian, which describes the low-energy excitations as 
excitations of semihuxons. The bifragmentation is readily explained in this limit, and 
the ground state of the system is found to be a macroscopic Schrodinger cat state of 
discrete semihuxons with opposite currents. 
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Appendix A. Diagonal structure of reduced density matrices 

The many-body ground state of a Bose-Hubbard Hamiltonian of N atoms in a triple-well 
potential, can be expressed as 

m = Y.Ck,i%i.N-k-i), 

k,l 

where the Fock basis of the system can be written as a product state of the reduced 
Fock basis {\k)j} for each subsystem j = 1, 2, 3: 

\k,l,N -k-l) = \k)i®\l)2®\N -k-l):^. 

The reduced density matrix of subsystem 1 can be computed as 


pi = '^2{m\3{n\ I ^ Ck,i\k,l,N - k - /) ^Q A - k' - l'\ j |m)2|n)3 

m,n \ k,l k',l' 

dk,m dN-k-l,n dN-k'-V,n ^k',m Ck',l'Ck,l 

m,n k,k' 1,1' 


\N — m — n){N — m — n\. 


(A.l) 


m,n 

From this expression one can see that the reduced density matrix of one subsystem i is 
diagonal in the reduced Fock basis of the corresponding mode {l/c)} = : 

N 


k=0 


(A.2) 


with {A^} the Schmidt coefficients that contain the information concerning the 
entanglement properties of subsystem i with the rest of the system. 
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